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A computation of raicleosynthesis bounds on the masses of long-lived Dirac and 
Majorana neutrinofil is reviewed. In particular an explicit treatment of the "differ- 
ential heating" of the Ve and Pe ensembles due to the residual out-of-equilibrium 
annihilations of decoupled heavy neutrinos is included. The affect is found to be 
considerably weaker than originally reported by Dolgov et al. p. For example, the 
bounds for a Dirac tau neutrino are rrii/^ < 0.37 MeV or TTLi/^ > 2S| MeV (for 
ANi/ > 1), whereas the present laboratory bound is mi/^ < 23.1 MeVcl. 



1 Introduction 

Nucleosynthesis considerations have proved to be an effective tool in finding 
limits on particle properties such as masses, couplings, lifetimes, neutrino mix- 
ing parameters and so on. The constraining power of nucleosynthesis sensi- 
tively depends on the constraints on primordial light element abundances that 
must be inferred from the observational evidence. This is the difficult part 
of NS considerations (see K.A. Olive, these proceedings). Nevertheless, as far 
as NS bounds on new physics are concerned, these details can be summarized 
by a single parameter: the number of equivalent neutrino degrees of freedom 
AiVy. The theoretical prediction for the helium abundance on the other hand, 
with or without new physics, is relatively straightforward and can usually be 
done very accurately. Parametrizing the deviation of the prediction from the 
standard result in units of AN^, one obtains a mapping of the nucleosynthesis 
bound on the space of parameters of the model. This is the situation in par- 
ticular for long-lived massive neutrinos, and this is what I mean with "precise' 
nucleosynthesis bounds. 

Massive annihilating particles affect the helium abundance indirectly by 
altering the expansion rate, or by directly altering the rate of reactions holding 
neutrons and protons in equilibrium: 

n + i^e ^ p + n + e'^^p + Pe, (1) 

during the time ~ 0.7 MeV when the n/p ratio is freezing out. Also, the 
changes in the expansion rate alter the time available for the free neutron decay 
n p + e + Je- 

For example, additional mass density speeds up the expansion rate relative 
to reactions (Q), which causes the n/p ratio to freeze out earlier, leaving behind 
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more neutrons and hence eventually more helium than in the reference case of 3 
massless neutrinos. Secondly, since 0{few) MeV neutrinos freeze out at rather 
low temperatures, their annihilations to Ve^e final states at temperatures below 
the chemical freeze-out temperature of v^s (Tchcm — 2.3 MeV) can produce 
some exess in the and Vf, number densities ("bulk heating"), leading to an 
increase of the overall rate of eqs. (Q). As a result equilibrium is maintained 
longer, leading to less neutrons and eventually less helium being produced. In 
section 2 I will review a computation that accurately accounts for these two 
effects. 

Moreover, the residual annihilations of already decoupled heavy neutrinos 
below Tkin — 1 MeV, when i/g's fall from kinetic equilibrium, j-jsan cause a 
deviation from equilibrium in the tail of the and Ve spectraa. Although 
small in amplitude, the effect of these neutrinos is boosted by the quadratic 
dependence on energy of reactions (^; the detailed balance argument relating 
the equilibrium conversion rates n p and p ^ n does not apply, and simply 
because there are more protons than neutrons around, the presence of this 
distortion has the tendency of increasing the relative number of neutrons and, 
eventually, the final helium abundance. I will compute the contribution to 
ATVi^ from this "differential heating" in section 3. All effects will be included 
in the final results, to be presented in section 4. 

2 Elements of the Computation 

The relevant momentum-dependent Bolzmann equations for the scalar phase- 
space distribution functions have the form: 

E,{dt - pHdp)Mp, t) = Ce,,(p, t) + t), (2) 

where Ei = [p^ + mfy^^ and H = (87rp/3M|[)^/^ is the Hubble expansion 
rate (with Mpi the Planck mass and p the total energy density). The index i 
runs over all particle species in the plasma; each distribution function of each 
species has an equation like (^) , and all of them are coupled through the elastic 
and inelastic collision terms C'e (p, t) and Ci {p, t) . 

Since the elastic scattering rates are much higher than the annihilation 
rates, neutrinos freeze out from the chemical equilibrium remaining in good 
kinetic contact with the other particles in the plasma. Their spectra can then 
be described to a very good accuracy by the pseudo-chemical potentials Zi(T), 
i.e. {Ei = yjp^ + mf) 

fip,z.) EE (e^'/^-+- + l)-i; T. ^ (^ 4 + 2Mr,) y/^ 
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Figure 1: la) Electron and tau neutrino number densities with a Majorana Ut with m^^ = 5 
MeV, normalized to equilibrium density no = {3({3) /2-k^)Tj^ . The dashed line for compari- 
son shows the unperturbed electron neutrino temperature, lb) Energy densities of Pt—, i't+ 
and Ue with a Dirac Ut with rriiy^ = 20 MeV, normalized to po = 77r^T^'*/240. Also shown 
is the energy density corresponding to Ut in the helicity averaged approximation. 



where he{T) is related to the entropy stored in the electrons and positrons: 
Se = 27r^/ie7^^/45. With the introduction of pseudo-chemical potentials the 
Bolzmann equations (H) can be integrated over the momenta. One then has a 
set of coupled ordinary differential equations for the pseudo-chemical potentials 
of each neutrino species and the photon temperature T-fJ. 

In fig. la a particular solution of the equation network with a mass to^^ — 5 
MeV is shown. The importance of tracking the Ve density is clearly visible: 
the number density of ly^s remains close to the equilibrium value until T ~ 2.3 
MeV, after which it freezes from chemical equilibrium. Since annihilations of 
i/r 's are still occurring at T ^ 2 MeV, a slight heating of the electron neutrino 
distributions results. In the units AN^, the corresponding effect ifl SAN^, ~ 
—3.6dn^^. In this example Srii,^ ~ 0.1 so that one expects to get SAN^, ^ —0.4, 
in good agreement with the full nucleosynthesis computatioiEJ. 

In the Dirac case there is the additional complication of different helicity 
states having different interaction strengths at high energies: positive helicity 
states will freeze out earlier with larger relic density (and the negative helicity 
states later with smaller) than the hypotetical Dirac neutrino with helicity- 
averaged interaction strength. One sees this effect in fig. lb, which shows the 
energy density stored into different helicities. These deviations almost exactly 
cancel each other in the sum, however, so that the total number density, and 
hence the induced effect on helium synthesis, is very close to the results found 
using the helicity-averaged approach. Figure lb also shows how the relative 
energy density stored into tau neutrinos sharply increases at small T. 
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3 Differential Heating 

Qualitatively thejahysics was explained in the introduction; for more details 
see Dolgov et al. Since one is considering a very small amplitude effect in 
the electron neutrino distributions, one can to a good accuracy use the lin- 
earized Bolzmann equation for the deviation. Moreover, because the effect 
is dominated by large energies, it is an exellent approximation to replace the 
Fermi-Dirac distributions by Bolzmann distributions. Then noting that the 
l.h.s. of eq. annihilates the equilibrium distribution, and trading the vari- 
ables {t,p) for Xi, = ra^^/Ti, and y = p/T^, one gets the generic equation for 
the deviation: 

d 

Hx^- — 5f{x, y) = Cann(a;, y) - Cci{x, y)Sf{x, y), (4) 
dxi, 

where x = m^^ /T . The variables x and Xu are of course simply related; using 
Xn allowed succinctly including the effect of entropy release in eq. (^. The 
annihilation term Cann(a^,y) acts as a source and it is peaked around y x. 
The term Cci{x,y) represents the restoring force of the elastic scatterings. In 
magnitude, Cd ^ Cann- Annihilation terma^are easily computed, and I find 
the forms similar to the ones by Dolgov et al.l3. On the Dirac case, for example 

Ca^„„ ~ 9.0 X 10-3^3^ ^^M!_^£aM! X cxp[-(x - y)yy]FD{x,y), (5) 

where n's are the number densities normalized as in fig. 1 and Foix^y) is a 
function with the asymptotic value of 1 for y ^ 0, a; ^ oo □. In numerical 
work I used exact functions, valid for any x and y, but this is not a significant 
effect. The elastic scattering term is given by (neq(O) = 1) 



The solution of is straightforward: 
^f{x,y)= / dx'A{x')Cann{x' ,y) exp 







dx"A{x")ai{x",y) 



(6) 
(7) 



where 



1 T-y / dhi 



and hi{T^) — 2 + he{T~,). Equations (^) and (Q) differ from those of Dolgov et 
al.a in that I included the dilution due to the entropy release in the electron 
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Figure 2: a) The distributions y'^&f[x, y) for m^^ = 10 MeV, and T = 1 MeV {x = 10) are 
shown with the corresponding bulk-equilibrium distributions, b) The effect on BBN of the 
differential heating in units of equivalent u degrees of freedom. 



annihilations, as well as the correct energy dependence in the expansion rate 
i?, both of which tend to weaken the effect. Moreover, the elastic scattering 
term (^) is somewhat larger than that of Dolgov et al., which also weakens the 
effect. Fig. 2a displays a particular solution of equation (Q) and final results 
for the effect of differential heating are shown in figure 2b. The effect is largest 
for the Dirac neutrino at around m^^ ~ 10 MeV. However, for m^^ — 25 MeV 
it is already below 0.25 effective neutrino degrees of freedom, i.e. a fourth of 
that found by Dolgov et al. These results are in very good agreement with the 
full numerical solution of the Bolzmann equation by Hannestad and MadserQ 
in the Majorana case. 

4 Results 

The mass bounds can be expressed in terms of fit functions of AiVi/ . Including 
all the effects to the electron neutrino distributions discussed above, we find 
that the Majorana masses are bounded by 

< v^(0.350 + 0.0477^+0.5912;) 61(0.15 -x) 
+ (0.079 + 0.576a; + 0.598a;2 - 0.514^^ + 0.211a;'*) e{x - 0.15) 
ml^ > 68.59- 63.83a; + 49.18a;2- 33.09x3 + 13.182;'*- 2. 17a;^ (9) 

where the units are MeV and I used x = ANi,. One sees that opening up a 
window for a stable tau neutrino below the laboratory bound of 23.1 MeVH 
would require relaxing the nucleosynthesis bound to AN^, > 1.44. Then the 
NS bound of ANi, < 1, together with the above laboratory bound rules out a 
long-lived Majorana tau neutrino with rn*j^ > 0.95 MeV. 
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In the Dirac case the upper bound on the disallowed region leads to the 
constraint 

m^^ > 41.88 - 28.47a; + 20.Ux^ - UMx^ + 5.49x'' - 0.90a;^ (10) 

where the units are again MeV. The effect of differential heating was for AiV^ > 
1 to increase the bound from 22 MeV to m^^ > 25 MeV, closing the window 
below the experimental bound of m^^ < 23.1 MeV at 95% CL. 

The computation of the lower bound on the disallowed region is entirely 
different from that in the preceding cases and completely unchanged by the 
differential heating effects. I will only quote the results from Fields et al.: for 
the bound AiV^, > 1 they give < 0.31 MeV and m^^ < 0.37 MeV, with 
Tqcd = 100 MeV. These constraints are conservative in the sense that they 
would-be somewhat stricter if the QCD transition temperature was chosen 
highefl. 

In conclusion, even with the weak constraint of AN^, < 1, nucleosynthesis 
bound is strong enough to exclude a long-lived (t > 100 sec) Majorana tau 
neutrino with m,y^ > 0.95 MeV and a Dirac tau neutrino with m^^ > 0.37 
MeV. 
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